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A summary of neutrino oscillation results is given along with a discussion of neutrino mass 
generation mechanisms, including high and low-scale seesaw, with and without supersym- 
metry, as well as recent attempts to understand flavor. I argue that if the origin of neutrino 
masses is intrinsically supersymmetric, it may lead to clear tests at the LHC. Finally, I 
briefly discuss thermal leptogenesis and dark matter. 

1 Status of neutrino oscillation experiments 

The discovery of neutrino oscillations provides the first evidence of physics beyond the Standard 
Model (SM), marking the beginning of a new era in particle physics. Thanks to their brilliant 
confirmation by reactor and accelerator experiments, oscillations constitute the only viable 
explanation for the observed flavor conversion of "celestial" neutrinos [TH3], requiring both 
neutrino mass and mixing, as expected in theories without conserved lepton number [4j[5]. 

Even in its simplest 3x3 unitary form, the lepton mixing matrix K = LO23UJ131-O12 differs 
from the quark mixing matrix in that each lu factor carries a physical phase: one is the KM- 
analogue and appears in oscillations, while the other two are Majorana phases and appear in 
lepton number (L)-violating processes. Current experiments are insensitive to CP violation, so 
that oscillations depend only on the three mixing angles ^12,^23,^13 and on the two squared- 
mass splittings Am'^i = m\ — ra\ and Am§]^ = m§ — characterizing solar and atmospheric 
transitions. Setting Am^x = in the analysis of atmospheric and accelerator data, and Am§]^ 
to infinity in the solar and reactor data analysis one obtains the neutrino oscillation parameters, 
as summarized in Figs. [1] and [21 Fig. [T] gives the allowed values of "atmospheric" and "solar" 
oscillation parameters, ^23 & Ajti^j^, and 612 & Am^x, respectively. The dot, star and diamond 
in the left panel of Fig. [T] indicate the best fit points of atmospheric, MINOS and global data, 
respectively. Similarly the "solar" oscillation parameters are obtained by combining solar and 
reactor neutrino data, as shown in the right panel. The dot, star and diamond indicate the 
best fit points of solar, KamLAND and global data, respectively. In both cases minimization 
is carried out with respect to the undisplayed parameters. One sees that data from artificial 
and natural neutrino sources are clearly complementary: reactor and accelerators give the best 
determination of squared-mass-splittings, while solar and atmospheric data mainly determine 
mixings. The right panel in Fig. [2] shows how data slightly prefer a nonzero ^13 value, though 
currently not significant, leading to an upper bound at 90%C.L. (3cr): 

0.060 (0.089) (solar+KamLAND) 
sin^ 6*13 < { 0.027 (0.058) (CH00Z+atm+K2K+MIN0S) (1) 
0.035 (0.056) (global data) 
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Figure 1: Current neutrino oscillation parameters, from a global analysis of the world's data [3] 
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Figure 2: Constraints on sin^ 6*13 from different neutrino oscillation data sets [3j. 



given for 1 dof, while the regions in Fig. [2] (left) correspond to 90% CL for 2 dof. The con- 
firmation of a non-zero ^13 would strongly encourage the search for CP violation in upcoming 

neutrino oscillation experiments [6l[7]. Note that the small parameter a = ^™"| is currently 

well-determined experimentally as a = 0.032 , 0.027 < a < 0.038 (3cr) . 

Before closing let us note that many effects may distort the "celestial" neutrino fluxes reach- 
ing our detectors, such as regular [5HTU] and random pUfT^ solar magnetic fields. These would 
induce spin-flavor precession in the convective zone |13H15j as well as density fluctuations deep 
inside the Sun's radiative zone [TBHTS] . Although these can modify the solar neutrino survival 
probabilities p"9ll20] , they can not have an important impact on the determination of oscillation 
parameters, thanks to the KamLAND reactor neutrino spectrum data. The result is that oscil- 
lations remain robust against astrophysical uncertainties, and of all oscillation solutions allowed 
by solar data [21], only the large mixing angle solution survives KamLAND's measurements [22] . 

Often the generation of neutrino mass in gauge theories (left panel in Fig. |4| is accompanied 
by effective sub- weak strength (^ eGp) flavour-changing (FC) or non-universal (NU) dimension- 
6 operators, as seen in the right panel. Such non-standard neutrino interactions (NSI) are 
expected in low-scale seesaw schemes, such as the inverse j23H25j and the linear [26] seesaw. In 
such schemes NSI would arise from the effectively non-unitary form of the corresponding lepton 
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mixing matrix [5] [65]. Relatively sizeable NSI strengths may also be induced in models with 
radiativcly induced neutrino masses [66p67j . Current determination of solar neutrino parameters 
is not yet fully robust against the presence of large NSI, allowing for a new "dark side" solution 
that survives the inclusion of reactor data [57] . 

In contrast, thanks to the large statistics of atmospheric data over a wide energy range, 
the determination of atmospheric parameters Amgj^ and sin^ 623 is fairly robust even in the 
presence of NSI, at least within the 2-neutrino approximation |28| . a situation likely to improve 
with future neutrino factories |29| . However, NSI operators may have dramatic consequences 
for the sensitivity to ^13 at a neutrino factory [30] and may affect the interpretation of future 
supernova neutrino data in an important way |31H33j . Improved NSI tests will also shed light 
on the origin of neutrino mass, helping discriminate between high and low-scale schemes. 

2 Neutrino mass and neutrinoless double beta decay 

Neutrino oscillations can not probe absolute neutrino masses, for this we need cosmic microwave 
background and large scale structure observations [34], high sensitivity beta decay and 
(neutrinoless double beta decay) studies [35] . The observation of neutrino oscillations suggests 
that light Majorana neutrino exchange will induce 0;//3/3 as illustrated in the left panel of 
Fig. [21 This nuclear process would hold the key to probe the nature (Dirac versus Majorana) 
of neutrinos |36| since, in a gauge theory, it would imply a Majorana mass for at least one 
neutrino [36], as illustrated in the middle panel of Fig. |3l 




Figure 3: Neutrino mass mechanism (left), black box theorem (center) [36] and OiyfUf] decay 
parameter versus lightest neutrino mass for inverse and linear seesaw models (right), from |37j . 

Such "black-box" theorem [35] holds in any "natural" gauge theory, though its implications 
are rather difficult to quantify in general [5S]- The 01/(3/3 detection prospects were discussed 
in |39j and are summarized in the right panel in Fig. [S] The broad bands are allowed by 
current neutrino oscillation data [2ll3] for normal and inverse neutrino mass hierarchy, while 
the narrow bands correspond to the case of lepton mixing angles fixed to the tri-bimaximal 
values [37]. Note that for normal hierarchy there is in general no lower bound on 01/(3/3 as there 
can be a destructive interference among the three neutrinos^. In contrast, the inverted neutrino 
mass hierarchy always gives a generic "lower" bound for the Ov/3/3 amplitude. On the other 
hand, quasi-degenerate neutrino models [15JI31] give the largest possible 0v/3(3 signal. Taking 
into account state-of-the-art nuclear matrix elements |45[ one can determine the best current 
limit, which comes from the Heidelberg-Moscow experiment, as well as future experimental 
sensitivities [33], summarized in the right panel in Fig. [5] for GERDA, Majorana and CUORE. 

^Specific flavor models may, however, lead to a lower bound on the Oi/fSf5 decay rate even for normal hierarchy 
neutrino spectra, as discussed in Refs. [37] I4UH42I . 
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3 Generating neutrino masses 

Underpinning the origin of neutrino mass remains a challenge despite the tremendous progress 
we have achieved. Neutrino masses arc markedly different from those of charged fermions in 
the SM. The latter arise by coupling the two chiral species to the Higgs scalar, hence linear 
in the electroweak symmetry breaking vacuum expectation value (vev) ($) of the Higgs scalar 
doublet ^ = H. By contrast, being electrically neutral, neutrinos may get mass with just one 
chiral species: in other words, on general grounds they are of Majorana type [S]. The lowest- 
dimensional lepton number violating (LNV) operator has AL = 2, namely AL3>L$, where L 
denotes a lepton doublet [i] , see left panel in Fig. HiEI- 




Figure 4: Neutrino mass [4] and non-standard neutrino interactions (NSI) operators [5]. 

Irrespective of their specific origin, the smallness of neutrino masses would come from the 
fact that they violate lepton number. The big issue is to identify which mechanism gives rise to 
the L-violating operator, its associated mass scale and its flavor structure. As for its magnitude, 
it may naturally be suppressed either by a high-scale Mx in its denominator, or may involve a 
low-mass-scale in its numerator. 

It is often argued that gravity breaks global symmetries (48l|49]. This would induce the L- 
violating operators suppressed by the Planck scale. The resulting Majorana neutrino masses are 
too small, hence the need for physics beyond the SM [SO] . It is usual to assume that this physics 
lies at a large sub-Planck scale, say, associated to unification. However, A may be suppressed by 
small scales, Yukawas and/or loop- factors [51]. There are three classes of mechanisms: (i) tree 
level, (ii) radiative, and (iii) hybrid, all of which may have high- or low-scale realizations, sug- 
gesting a fair chance that the origin of neutrino mass may be probed at accelerator experiments 
like the LHC. Depending on the nature of spontaneous lepton- number breaking there may be 
an extra neutral gauge boson HSJIH] or a Nambu-Goldstone boson coupled to neutrinos [S3]. 

3.1 Minimal high-scale seesaw 

Weinberg's dimension-5 operator [J] may arise from the exchange of heavy fermion states with 
masses close to the "unification" scale. Depending on whether these are SU{2) singlets or 
triplets the mechanism is called type-I [5^H57] . or type-Ill seesaw [SS], respectively. Neu- 
trino masses may also arise from the exchange of heavy triplet scalars, now called type-II 
seesaw [5] j53l[59| , as seen in the right panel in Fig. [5l The hierarchy of vevs required to account 
for the small neutrino masses <^ V2 ^ vi is consistent with the minimization of the scalar 
potential. The resulting perturbative diagonalization of the seesaw mass matrix was given in 
Ref. |53j in the most general form that may be used in any model. From a phenomenological 
viewpoint, however, the most basic effective description of any seesaw is in terms of the SM 
SU{3) (g) SU{2) (g) U{1) gauge structure with cxphcit L- violation [Hj. 

^Note that neutrino masses may also arise from higher dimension effective operators I46II47I . 
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Figure 5: Type-I and III (left) and Type-II (right) realizations of the seesaw mechanism. 

3.2 "Non-minimal" seesaw 

The seesaw mechanism is not any particular theory but rather represents a broad language in 
terms of which to phrase neutrino mass model-building and many are its pathways. Indeed, 
it can be implemented non-minimally with lepton number broken explicitly or spontaneously, 
over a wide range of energy scales, in a variety of models with different gauge groups and 
multiplet contents, with or without supersymmetry. There is no point to attempt giving here 
a full taxonomy of seesaw schemes. However one should stress that any model must ultimately 
reduce to the SM. Hence what is phenomcnologically most relevant, for example to describe 
neutrino oscillation data, is the effective structure of the SM seesaw mixing matrix, given in [5] . 
In addition to the mixing angles characterizing oscillations, the latter includes non-unitarity 
effects which will become an important topic in the agenda of future studies probing neutrino 
propagation beyond oscillations. 

An attractive class of non-minimal seesaw schemes employs, in addition to the left-handed 
SM neutrinos vl, two SU{3) «) SU{2) U{1) singlets v", S [23] (see also e.g. [6lH64|). The 
basic parameter characterizing the violation of lepton-number can be small [25 1 137 ] and may be 
calculable from supersymmetric renormalization group evolution effects [24j . 

3.3 Radiative neutrino masses 

Neutrino masses may be radiatively calculable [SBJIHI], with no need for a large scale. In this 
case the suppression comes from small loop-factors and Yukawa couplings. The same way as 
low-scale seesaw schemes, the states responsible for generating neutrino masses in radiative 
models may lie at the weak-scale, opening the door to phenomenology at the LHC [60j . 

4 Understanding and probing flavor 

There is no reasonable doubt that flavor is violated in neutrino propagation [2l[3|. Current 
oscillation experiments indicate solar and atmospheric mixing angles which are unexpectedly 
large when compared to quark mixing angles. To a good approximation they are given by [68], 

tan^ 6*23 = tan^ 6^3 = 1 sin^ 6'chooz = sin^ e'Jg = tan^ ^sol = tan^ 9^^ 0-5- (2) 

Understanding the pattern of lepton mixing angles from first principles constitutes a big 
challenge to unified theories of flavor where quarks and leptons are related. Many less ambitious 
schemes have been suggested in order to reproduce the tri-bi-maximal pattern, at least partially, 
using various discrete flavor symmetry groups containing mu-tau symmetry, e. g. [43114411691176] . 
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In general one expects the flavor symmetry to be valid at high energy scales. Deviations from 
tri-bi- maximal ansatz |77j may be calculable by rcnormalization group evolution |78H80j . A 
simple possibility is that, as a result of a given flavor symmetry such as A4 [43l|44], neutrino 
masses unify at high energies [5T], the same way as gauge couplings do, due to supersymmetry. 
Such quasi-degenerate neutrino scheme predicts maximal atmospheric angle and vanishing ^13, 
^23 = 7r/4 and ~ , leaving the solar angle 612 unpredicted, but Cabibbo-unsuppressed, 
O12 = If CP is violated ^13 becomes arbitrary and the Dirac phase is maximal [70]. The 

lower bound on the absolute Majorana neutrino mass scale mp > 0.3 eV ensures that the model 
will be probed by future cosmological data and 0;//3/3 searches. 

It is natural to expect that, at some level, lepton flavor violation will also show up as 
transitions involving the charged leptons, since these sit in the same electrowcak doublets as 
neutrinos. Rates for lepton flavour violating processes Ij k + j often lie in the range of 
sensitivity of coming experiments, providing an independent test. There are two basic mecha- 
nisms: (i) neutral heavy lepton exchange and (ii) supersymmetry |5SI - [57| . Both exist 
in supersymmetric seesaw- type schemes of neutrino mass, the interplay of both depends on the 
seesaw scale, and has been considered in [88]. Barring fine-tunings, high-scale seesaw models 
require supersymmetry in order to have sizeable LFV rates. The most interesting feature of 
these models is that they bring in the possibility of direct lepton flavour violation in the pro- 
duction of supersymmetric particles. As seen in Fig. |5] this provides the most direct way to 
probe LFV at the LHC in high-scale seesaw models. 




[GeV] [GeV] 

Figure 6: LFV rate for fi-r lepton pair production from X2 decays versus M1/2 for the indicated 
Too values, assuming minimal supcrgravity parameters: fi> 0, tan/3 = 10 and Aq = GeV, for 
type-I (left) and for type-II seesaw (right). Here Ai — 0.02 and A2 = 0.5 are Type-II seesaw 
parameters, and we imposed the contraint Br(/^ ^ e + 7) < 1.2 ■ 10^^^, from Ref. [89]. 

In low-scale seesaw schemes, by contrast, the sizeable admixture of "right-handed" (RH) 
neutrinos in the charged current (rectangular nature of the lepton mixing matrix [5]) induces 
potentially large LFV rates even in the absence of supersymmetry |82j . Indeed, an important 
point to stress is that LFV [82l|83] and CP violation [9T1|92| can occur in the masslcss neutrino 
limit, hence their attainable magnitude is unrestricted by the smallness of neutrino masses. In 
Fig. [7] we display Br{ii e-f) versus the small LNV parameters /i and for two different 
low-scale seesaw models, the inverse and the linear seesaw, respectively. Clearly the LFV rates 
are sizeable in both cases, the different slopes with respect to /i and vl follow from the fact 
that AL = 2 in the first case and AL = 1, in the second. 



6 



LP09 




Figure 7: Left: Br{ii — > ej) versus the LNV scale for inverse seesaw (top: red color) and linear 
seesaw (bottom, blue color). In both low-scale seesaw models we fix M = 100 GeV^ (continuous 
line), M = 200 GeF (dashed line) and M = 1000 GeV^ (dot-dashed line), from [37]. The right 
panel shows typical correlation between mu-e conversion in nuclei and Br{p — > 67), from [90] . 



Similarly |90j in low-scale seesaw models the nuclear /x~ — e~ conversion rates lie within 
planned sensitivities of future experiments such as PRISM [S3] . Note that models with specific 
flavor symmetries, such as those in [IHIISZ] relate different LFV rates. To conclude we mention 
that some seesaw schemes, like type-Ill with tiny Yukawas [HH] or inverse type-Ill [35], may be 
directly probed at the LHC by directly producing the TeV states with gauge strength [60) . 

5 Probing neutrino properties at the LHC 

In supersymmetric models lepton number can be broken together with the so-called R parity, 
leading to an intrinsically supersymmetric origin for neutrino masses |94H96| . This may happen 
spontaneously, driven by a nonzero vev of an SU{3) (g) SU (2) (g) U (1) singlet sneutrino |S7l - [^ . 
leading to an efi^ective model with bilinear violation of R parity |1001ll01j . The latter provides 
the minimal way to break R parity and add neutrino masses to the MSSM (lOlj . One finds 




200 ■100 eoo soo 1000 1200 1400 1600 1800 2000 



Figure 8: Xi decay length versus mo for Aq = —100 GeV, tan/3 ~ 10, fi > 0, and various mi/2 
values. The three shaded bands around mi/2 = 300, 500, 800 GeV correspond to the variation 
of the BRpV parameters in such a way that the neutrino masses and mixing angles fit the 
required values within 3cr. The right panel gives the ratio of branching ratios, Br(x5 — J* Z^?'?) 
over Br(x5 'tq'q) in terms of the atmospheric angle in bilinear R parity violation |103j . 

that, typically, the atmospheric scale is generated at tree level by neutralino-exchange weak- 
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scale seesaw, while the solar scale is radiatively induced |102j . Unprotected by any symmetry, 
the lightest supersymmetric particle (LSP) decays. Given the masses indicated by neutrino 
experiments these decays will happen inside typical detectors |102H104j but with a decay path 
that can be experimentally resolved, leading to a so-called displaced vertex jl05| , see left panel 
in Fig. |8l More strikingly, LSP decay properties correlate with the neutrino mixing angles. In- 
deed, as seen in the right panel in Fig. [5] the LSP decay pattern is predicted by the low-energy 
measurement of the atmospheric angle [1031ll06lll07j . allowing for a clear test at the LHC [50] . 
namely a high-energy redetermination of ^23- Similar correlations hold in variant models based 
on alternative supersymmetry breaking mechanisms, where other states appear as LSP [108j . 

6 Neutrinos as cosmological probes 

Neutrinos can probe very early epochs in the evolution of the Universe, previous to the elec- 
troweak phase transition. For example, the high-scale generation of neutrino masses through 
the seesaw mechanism may induce the observed baryon asymmetry of the Universe, as well as 
the dark matter, as I now discuss. 

6.1 Thermal leptogenesis 

The observed cosmological matter-antimatter asymmetry in the Universe may arise from the 
C/CP-violating out-of-cquilibrium decays of the heavy RH neutrinos present in the seesaw. 
These take place before the clcctrowcak phase transition |109j through the diagrams in the 
left panel in Fig. [S] The Icpton asymmetry thus produced gets converted, through sphaleron 
processes, into a baryon asymmetry. This so-called leptogenesis scenario [llOlllllj is a frame- 



Figure 9: Diagrams (left), flavor effects on minimum required leptogenesis scale (right) jll3j . 

work to explain just one number, namely the baryon asymmetry, currently well-determined by 
WMAP |112j . It turns out that, as displayed in the right panel in Fig. [HI from Ref. }113] . after 
taking into account carefully washout and flavor effects one finds that the required asymmetry 
can be achieved for a given range of model parameters which fits neutrino masses indicated by 
neutrino oscillation experiments. It did not have to be so, a priori, hence this may be taken as 
a success. The figure also shows how the inclusion of fiavor effects lowers the minimum value 
of the lightest RH neutrino mass required for successful leptogenesis. Nevertheless, in order 
to avoid gravitino overproduction jll4j . which would destroy the standard Big Bang Nucle- 
osynthesis predictions, one also requires an upper bound on the reheat temperature Tn after 
inflation, incompatible with Fig. [9] |115j . One way to prevent such gravitino crisis is to assume 
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enhancement coming from resonant leptogenesis |116) . Alternatively, there are many ways to 
go beyond the minimal type-I supersymmetric seesaw |117H119] . For example one may add a 
small R-parity violating superpotential term XiV^iHuHd, where v'^i are RH neutrino supermulti- 
plets |118| . In the presence of this term the produced asymmetry can be enhanced. In extended 
SO(IO) supersymmetric seesaw schemes leptogenesis can occur at relatively low scales, through 
the decay of a new singlet, as illustrated in the left panel in Fig. [TUl This not only avoids the 
gravitino crisis but also opens the possibility of detecting the new neutral gauge boson at the 
LHC |1191I120| . The right panel illustrates how a sizeable asymmetry may be achieved just 
with the leptonic CP violation parameter 5 that characterizes neutrino oscillations. 




Figure 10: Low-scale leptogenesis in supersymmetric SO(IO) models, from |119|ll20j . 



6.2 Neutrino masses and dark matter 

Neutrinos may get mass through the spontaneous breaking of ungauged lepton number. Due 
to quantum gravity effects the associated Goldstone boson - the majoron - is likely to pick up a 
mass, and play the role of late-decaying Dark Matter, decaying mainly to neutrinos |1211ll22j . 
Cosmic microwave background observations place constraints on the majoron lifetime and mass, 
illustrated in left and middle panels of Fig. 1111 This decaying dark matter scenario arises in 




Figure 11: Late decaying majoron dark matter: decay parameters allowd by the CMB |123j 
(left- middle); probing sub-leading decay to two photons J — > 77 (right), from |124j . 

type-II seesaw models, where the majoron couples to photons through the Higgs triplet and 
may be tested through the mono-energetic emission line from its sub-dominant decay to two 
photons, as illustrated in the right panel in Fig. [11] 

Neutrino masses may also open new possibilities for "conventional" supersymmetric dark 
matter. For example, within the inverse seesaw mechanism minimal supergravity is more likely 
to have a sneutrino as the lightest superparticle than the conventional neutralino. Such schemes 
naturally reconcile the small neutrino masses with the correct relic sneutrino dark matter abun- 
dance and accessible direct detection rates in nuclear recoil experiments [125) . 
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